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Abstract. We investigate carbon isotope effects on the infrared bands of a laboratory analogue of carbonaceous 
dust, the quenched carbonaceous composite (QCC), synthesized from a plasma gas of methane with various 
12 C/ 13 C ratios. Peak shifts to longer wavelengths due to the substitution of 12 C by 13 C are clearly observed in 
several absorption bands. The shifts are almost linearly proportional to the 13 C fraction. New features associated 
with 13 C are not seen, indicating that the infrared bands in the QCC are not very localized vibration modes 
but come from vibrations associated with rather large carbon structures. An appreciable peak shift (AA ~ 0.23- 
0.26 fim per 13 C fraction) is detected in the 6.2 /im band, which is attributed to a carbon-carbon vibration. A peak 
shift (AA ~ 0.16-0.18 /im per 13 C fraction) in an out-of-plane bending mode of aromatic C— H at 11.4 ^m is also 
observed, while only a small shift (AA < 0.015 /im per 13 C fraction) is detected in the 3.3 /im band, which arises 
from a C— H stretching mode. The present experiment suggests that peak shifts in the unidentified infrared (UIR) 
bands, particularly in the 6.2 /im band, should be detectable in celestial objects with low 12 C/ 13 C ratios (< 10). 
The isotopic shifts seen in the QCC are discussed in relation to the variations in the UIR band peaks observed 
in post-asymptotic giant branch stars and planetary nebulae. The observed peak shift pattern of the UIR bands 
is qualitatively in agreement with the isotopic shifts in the QCC except for the 7.7 /im band complex although 
the observed shifts in the UIR bands are larger than those inferred from derived isotope ratios for individual 
objects. The poor quantitative agreement may be attributed partly to large uncertainties in the derived 12 C/ 13 C, 
to possible spatial variations of the isotope abundance within the object, and to combinations of other effects, 
such as hetero-atom substitutions. The present investigation suggests that part of the observed variations in the 
UIR band peaks may come from the isotopic effects. 

Key words, infrared: ISM - dust, extinction - stars: AGB and post-AGB - planetary nebulae: general - ISM: 
abundances — ISM: lines and bands 



1. Introduction 

A wide range of the 12 C/ 13 C ratio has been reported in 
various celestial objects. The ratio changes due to nucle- 
osynthesis and mixing in the interior of stars. During the 
first dredge-up on the red giant branch (RGB) the con- 
vective envelope reaches regions abundant in 13 C that was 
processed from 12 C and the ratio decreases. Observations 
of RGB stars often show the 12 C/ 13 C ratio of 5-20, even 
lower than theoretically predicted, suggesting the presence 
of extra-mixing below the convective envelope (e.g. Gilroy 
1989}). In the third dredge-up during the asymptotic giant 
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branch (AGB) phase, an increase in 12 C/ 13 C is generally 
expected, but the ratio could also decrease due to cool 
bottom processing for low-mass stars (Wasserburg et at. 
IT§9"5l Nollett et al. l2Utj5) l or to hot bottom burning for 
more massive stars (Frost et al. D.998). Observations of 
five carbon-rich circumstellar envelopes indicate the ratio 
of 30-65 (Kahane et al. 1992), while ten carbon stars are 
shown to have the ratio in the range 12-60 in their circum- 
stellar envelopes (Greaves & Holland ll997f) . Some carbon 
stars show the ratio of 12 C/ 13 C as low as 3 (Lambert et 
al.lUSa Ohnaka & Tsufi lTMBl Schoier & Olofsson 123131) ■ 
The 12 C/ 13 C ratio in post-AGB stars and planetary neb- 
ulae (PNe) reflects the cumulative effects of different mix- 
ing and nuclear processing events during the entire evo- 



2 



Wada, Onaka, Yamamura et al.: Carbon Isotope Effects on IR Bands 



lution of their progenitors. Lower limits of the ratio of 3 
to 10 have been obtained for several objects in the post- 
AGB phase (Palla et al. EUUUl Greaves & Holland ITW7|l . 
Recently Josselin & Lebre lj2~0DTJl estimated an upper limit 
of 12 C/ 13 C of 5 for the post-AGB candidate, HD179821, 
whereas a relatively large ratio of 72±26 is reported for an- 
other post-AGB star, HD56126 (Bakker & Lambert ITM5|l . 
Clegg et al. (|1997[l found low ratios of 15 and 21 in two 
PNe. Further low values of the ratio of 2-30 have been re- 
ported in recent studies of several PNe (Palla et al. pOOO 
Balser et al. [2002 Josselin & Bachiller 2003), suggesting 
that some stars undergo non-standard processing in the 
stellar interior and a low 12 C/ 13 C can be expected during 
the late stage of their evolution. The solar system value is 
89 (Anders & Grevesse \TM% . 

The 12 C/ 13 C ratio also provides key information on 
the chemical evolution in the Galaxy (for a review, Wilson 
1999). Observations of molecules and solid CO2 in inter- 
stellar medium indicate that the ratio ranges from 10- 
100 and increases with the Galactocentric distance. The 
12 C/ 13 C ratio is suggested to be about 10-20 in the 
Galactic center region (Wilson 119991 Boogert et al. [2000 
Savage et al. 2002 ). The Galactic gradient is thought to be 
built by nucleosynthesis of the Galactic chemical evolution 
and the suggested ratio of 10-20 in the Galactic center in- 
dicates the presence of significant stellar sources of 13 C. 
Interstellar graphite spherules in the Murchison meteorite 
show a range of the ratio of 7-1330 (Bernatowicz et al. 
1999|. Some presolar SiC grains show very low 12 C/ 13 C 
ratios of less than 10 and they are thought to originate 
from very 13 C-rich stars in the AGB phase (Amari et al. 
2001). These observations suggest that low 12 C/ 13 C ratio 
environments are not uncommon in objects in the AGB, 
post-AGB, and PN phases as well as in some interstel- 
lar medium. Carbon-bearing species formed in these envi- 
ronments could thus show non-negligible carbon isotopic 
effects in their spectrum. 

A set of emission bands at 3.3, 6.2, 7.6-7.8, 8.6, and 
11.2 /im have been observed in various celestial objects and 
are called the unidentified infrared (UIR) bands. Fainter 
companion bands are also sometimes seen. The exact na- 
ture of the carriers has not yet been understood com- 
pletely, but it is generally believed that the emitters or 
emitting atomic groups containing polycyclic aromatic hy- 
drocarbons (PAH) or PAH-like atomic groups of carbona- 
ceous materials, including such as nanodiamond grains, 
are responsible for the UIR bands (Leger & Puget 119841 
Allamandola et al. 119851 Sakata et al. 1984 Papoular et 
al. rn?g§l Arnoult et al. lMTUl Jones & d'Hendecourt l2UUU|) . 
Alternatively, Holmid (2000) has recently proposed de- 
excitation of Rydberg matters as possible carriers. The 
UIR bands have been observed in a wide range of objects, 
including H II regions, reflection nebulae, post-AGB stars, 
and PNe (for a review, see Tokunaea I1997f) . They have 
also been commonly seen in the diffuse Galactic emission 
(Tanaka et al. lT§9l)l Onaka et al U^M Mattila et al. lT§9l)l 
Kahanpaa et al. 120031 as well as in external galaxies (e.g. 
Mattila et al. nMfl Helou et al. OHTTni Reach et al.^OOO 



Lu et al. I2003f) . indicating that the carriers are a com- 
mon member of interstellar medium and present in var- 
ious environments. Carbon-rich objects in the evolution- 
ary stage from post-AGBs to PNe often show the emis- 
sion bands and thus isotopic effects should be detectable 
if they arise from carbonaceous materials of low 12 C/ 13 C 
ratios. Simple calculations of a 13 C-benzene molecule sug- 
gest that the peak shift can be as much as 0.15 jam for the 
C— C stretching mode (AppendixEJ. 

Observations of the Infrared Space Observatory (ISO; 
Kessler et al. I1996fl have provided a large database of the 
UIR band spectra in various objects (e.g. Beintema et al. 

Molster et al. IT59l)l Verstraete et al. IT§9l l2DtjTl 
BoulangerEHHEl Cesarsky et al. I2000al l2000bl Uchida et 
al. IT5M1 EHCTil Moutou et al. HUM Honey et al. HE]}. 
Recently Peeters et al. (12002) have investigated in de- 
tail the 6-9 /Ltm spectra of 57 sources taken by the Short 
Wavelength Spectrometer (SWS; de Graauw et al. I1996f) 
on board the ISO and found that the 6.2, 7.7, and 8.6 
/im UIR bands show appreciable variations particularly 
for post-AGB stars and PNe. On the other hand, the vari- 
ations in the 11.2 /im band are relatively modest and those 
in the 3.3 /mi are less pronounced (Tokunaga et al. 119911 
Roche et al. 119911 Hony et al. 120011 van Diedenhoven et al. 
2003|. These variations can be interpreted in part by ni- 
trogen substitutions in PAHs and anharmonicity, but not 
all of the observed aspects of the UIR bands have yet been 
fully understood (Verstraete et al. 120011 Pech et al. lMEZl 
Peeters et al. 20Q2J . Part of the observed variations could 
also originate from isotopic effects of the UIR band carri- 
ers since the objects that show the variations are mostly 
post-AGB stars and PNe, in which small 12 C/ 13 C ratios 
can be expected. 

In the present paper we investigate isotopic effects on 
the UIR bands experimentally. We synthesize a labora- 
tory analogue of carbonaceous dust, the quenched car- 
bonaceous composite (QCC; Sakata et al. 1984), with var- 
ious 12 C/ 13 C ratios from the starting gas of a mixture of 
12 CH 4 and 13 CH4. The QCC shows infrared bands simi- 
lar to the UIR bands and shifts in the band peaks due to 
the 13 C substitution are clearly detected. In Sect. 2 we de- 
scribe the experimental procedure. The results are shown 
in Sect. 3 and discussed in comparison with observations 
in Sect. 4. A summary is given in Sect. 5. 

2. Experimental 

The experimental procedure for synthesizing the QCC 
is described in detail in Sakata et al. Q1984[l . Methane 
(CH4), the source gas of the QCC, is decomposed by 
the imposed microwave radiation and becomes a plasma. 
Carbonaceous condensates are formed in the injection 
beam of the plasma by quenching of the gas. Typically 
two types of the QCC are formed. One is a brown-black 
material (hereafter called dark-CC), which is collected on 
a substrate in the main injection beam. It has been shown 
to consist of a coagulation of carbon-onion-like particles 
(Wada et al. I1999J) . The other is a yellow-brownish mate- 
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rial (hereafter filmy-QCC), deposited on the surrounding 
region of the injected plasma gas (Sakata et al. I1987f) . It 
is a material rich in organic molecules, such as PAHs. We 
prepare the starting gas of CH4 with the 13 C fraction of 
1%, 11%, 25%, 45%, 65% and 99% by mixing 99% 13 C 
methane with the natural isotope abundance (1% 13 C) 
methane. We assume that the 13 C fraction in the QCC 
is equal to that of the starting gas in the present analy- 
sis because the reactions in the present experiments take 
place at high temperatures. 

The QCC is collected either on a KBr (for filmy-QCC) 
or on a BaF2 crystal (for dark-QCC) and the absorption 
spectra are measured at room temperature. The spec- 
tra are taken by a Perkin-Elmer 2000-FTIR spectrometer 
with the resolution of 4 cm" 1 . The spectra of the dark- 
QCC are measured after washing it with a small amount of 
acetone and removing organic molecules in the dark-QCC 
in a similar manner to previous experiments (Sakata et al. 

3. Results 

The isotopic substitution gives rise to a shift in the wave- 
length of vibration modes due to mass effects. With the 
replacement of 12 C with 13 C, the infrared band is expected 
to be shifted to longer wavelengths. 

3.1. Spectra of filmy-QCC 

The filmy-QCC is a coagulation of many kinds of organic 
molecules. Mass analysis of the gases evaporating from 
the filmy-QCC shows that they contain various kinds of 
PAHs. The filmy-QCC also emits red fluorescence under 
ultraviolet irradiation (Sakata et al. I1992fl . 

The spectra of the filmy-QCC with various 13 C frac- 
tions are shown in Fig. ^ They show several absorption 
features similar to the UIR bands observed in celestial ob- 
jects (Sakata et al. CHHll EHHII USHDJ - For the filmy-QCC 
with 1% 13 C the bands of the C— H stretching character- 
istic group appear at 3.29, 3.42, 3.51, and 3.53 /im (see 
also Fig.|5J). An aromatic C=C bond vibration appears at 
6.2 /im, and aromatic C— H out-of-plane bending modes 
are observed at 11.4, 11.9, and 13.2 /im. A weak feature 
also appears clearly at 7.6 /im in addition to the bands at 
7.0 and 7.3 /im. A very broad band is seen around 8.6 /im, 
but its peak wavelength cannot be determined accurately 
because of the weak and broad feature. All the features 
observed in the samples of various 13 C fractions have cor- 
responding features in the 1% 13 C sample and the peak 
wavelengths are shifted to longer wavelengths with the 13 C 
fraction. No new absorption band associated with 13 C ap- 
pears. 

A large shift (AA - 0.23 /mi between 1% to 99% 13 C) 
is observed in the 6.2 /im band. The C— C mode involves 
at least two carbon atoms and makes the shift large. The 
shifts in the peaks arising from C— H out-of-plane bending 
modes are also seen. The 11.39 /im peak at 1% 13 C shifts 
to 11.55 /im at 99% 13 C. On the other hand, the shift 
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Fig. 1. Absorption spectra of filmy-QCC with various 13 C 
fractions. The 13 C fractions are (a) 1%, (b) 11%, (c) 25%, 
(d) 65%, and (e) 99%. Each curve is shifted arbitrarily in 
the vertical direction. 
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Fig. 2. Spectra of the filmy-QCC in the 3 /mi region. The 
13 C fractions are (a) 1%, (b) 11%, (c) 25%, (d) 65%, and 
(e) 99%. Each curve is shifted arbitrarily in the vertical 
direction. 

in the 3.29 /xm band of aromatic C— H stretching is quite 
small (AA ~ 0.014 /xm). Fig.[2]shows the spectra expanded 
in the 3 /im region. It is noticeable that the profile of the 
3.3 /im band becomes narrower with the 13 C fraction. The 
cause of this change is unknown at present. A small shift 
is also observed in the 3.4 /im band, which is attributed to 
an asymmetric vibration of methylene (C— H2). 

The peak wavelength of the small feature at 7.6 /im 
is also shifted to longer wavelengths with the 13 C frac- 
tion (AA ~ 0.22 /im). The large shift indicates that the 
vibration involves more than one carbon atoms, such as 
C— C. The broad bump at 8.6 /im does not shift clearly. 
Sakata et al. (1987) assigned the 7.6 and 8.6/im peaks of 
the filmy-QCC to vibration modes of a kind of the ketone 
bond C=C— C=0. The 7.6 /im peak is always observed in 
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Fig. 3. Absorption spectra of (a) dark-QCC and (b) filmy- 
QCC with the natural isotopic fraction ( 13 C = 1.1%). 



the filmy-QCC spectrum although the relative strength 
changes slightly with individual runs. The exposure to air 
does not increase the strength of the 7.6 (im band. The 
present experiment thus suggests that the 7.6 fim compo- 
nent is a molecular product of the plasma gas and is not 
directly related to oxidation. 

3.2. Spectra of dark-QCC 

High-resolution electron-microscopy reveals that the dark- 
QCC is comprised of a coagulation of onion-like parti- 
cles of the diameter of 10-15 nm (Wada et al. 11999(1 . 
It is insoluble in acetone. Even after washing with ace- 
tone, the dark-QCC still contains some organic molecules. 
Mass spectroscopy indicates that they are mainly com- 
posed of compact PAHs. The dark-QCC adsorbs organic 
molecules that are formed together and also contains hy- 
drogen atoms in its structure. C— H bonds in both com- 
ponents will give rise to vibration modes in the infrared. 
Fig-Elshows the absorption spectra of the dark-QCC with 
1% 13 C fraction together with the filmy-QCC spectrum 
for comparison. The dark-QCC shows a strong continuum 
component, which decreases with the wavelength. Small 
peaks are observed at 3.3, 3.42, 5.8, 6.3, 7.0, 7.3, 8.2, 11.4, 
12.0, and 13.2 /mi. The 3.3, 3.42, 11.4, 12.0, and 13.2 /mi 
bands are attributed to CH vibrations, and the 6.3 /mi 
band is ascribed to an aromatic C=C bond vibration. The 
broad 8.2 fim band becomes strong when it is exposed to 
air. An increase of the 5.8 fim band is also seen during the 
exposure to air, suggesting that both of them are associ- 
ated with C— O— and C=0 bonds formed by oxidation. 

Fig. 01 shows the absorption spectra of the dark-QCC 
with various 13 C fractions. The continua were fitted with 
cubic spline functions and have been subtracted to show 
the absorption features clearly. As in the filmy-QCC spec- 
tra, peak shifts with the 13 C fraction are clearly observed 
and new features without corresponding ones in the 1% 
13 C sample are not seen in large 13 C fraction samples. 
A large shift (AA ~ 0.26 /mi) is again observed for the 
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Fig. 4. Absorption spectra of the dark-QCC with various 
13 C fractions with the continuum subtracted. The 13 C 
fractions are 1% (a), 11% (b), 25 % (c), 65% (d), and 
99% (e). Each curve is shifted arbitrarily in the vertical 
direction. 
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Fig. 5. Spectra of the dark-QCC in the 3 /mi region. The 
13 C fractions are (a) 1%, (b) 11%, (c) 25%, (d) 65%, and 
(e) 99%. Each curve is shifted arbitrarily in the vertical 
direction. 



6.3 //m band. The peak position of the 8.2 /mi band can- 
not be measured accurately because of the wide profile, 
but the peak is clearly shifted to longer wavelengths with 
the 13 C fraction. The band width also seems to change. 
Particularly for the 99% 13 C sample the 8.2 /im feature 
becomes very broad. The 3.3 /im band is quite weak in 
the dark-QCC and the peak position of the 3.3 /j,m is dif- 
ficult to determine (see Fig. 03). The feature almost fades 
away in the 99% 13 C sample. The change in the profile of 
the 3.3 /xm seen in the filmy-QCC sample is not observed 
in the dark-QCC. 



Wada, Onaka, Yamamura et al.: Carbon Isotope Effects on IR Bands 5 




ZO AO 60 30 100 ZO 40 60 30 100 

Fig. 6. Peak wavelength versus 13 C fraction of the filmy-QCC (open circles) and the dark-QCC (filled circles), (a) 
Aromatic C— H stretching modes, (b) aromatic C=C vibration modes, (c) 7.6 and 8.2 (im bands, and (d) C— H out- 
of-planc bending modes. Least-squares fit lines are also plotted. 

Table 1. Peak shift of the absorption bands in the QCC and the UIR bands 



Peak wavelength 


Slope of the fitted line 


Estimated shift between 


Observed range of 


of 12 C QCC 


(jim per 


13 C fraction) 


12 C and 13 C benzene" 


the peak wavelength 


(fan) 


filmy-QCC 


dark-QCC 


(fim) 


(pm) 


3.3 


0.014 ±0.001 


0.0063 ± 0.0007 c 


0.013 


3.288 - 3.297 


6.2 


0.231 ± 0.018 


0.262 ±0.023 


0.14 


6.20 - 6.27 (6.30) 
7.72 - 7.97 d 


7.6 


0.219 ±0.010 






8.2 




0.379 ±0.124 c 




(~ 8.22) 


11.4 


0.163 ±0.007 


0.183 ±0.007 


0.04 


11.20 - 11.25 



a See Appendix 1X1 

b The observed range between class A and B objects taken from Tokunaga et al. (199D, Peeters et al. 120021k and van 
Diedenhoven et al. (2003). The peak wavelengths of class C objects are indicated in parentheses (see text), 
c The slopes have large uncertainties because of the broad feature. 
d The peak wavelength range of the 7.8 fan subcomponent (see text). 



3.3. Peak shifts of the absorption bands 

Fig. El plots the peak wavelength of each band against 
the 13 C fraction. The change of the peak wavelength is 
roughly linear with the 13 C fraction for all the bands plot- 
ted. Linear lines obtained by least-squares fits are also 
shown. The slope of the fitted line is given in Tabled to- 
gether with the shifts estimated from the calculation of 
12 C- and 13 C-benzene (Appendix^). The amount of the 
shift in the 6.2 /im band is large and even larger than 
the estimate based on the simple calculation of benzene 
molecules. The shift of the 11.4 fim band is also larger than 
the calculation. The shift in the 3.3 fim band of the filmy- 



QCC is in agreement with the calculation, while that of 
the dark-QCC is smaller than the calculation, ft should 
be noted, however, that the 3.3 ^m band of dark-QCC is 
weak and the peak position cannot be determined accu- 
rately (Fig. [5J • 

The present experiment shows no new features cor- 
responding to 13 C. If the vibration mode is local, there 
should appear peaks corresponding to 12 C and 13 C. The 
spectral resolution is sufficiently high except for the 3.3 fim 
band, which shows only small shifts, but any appreciable 
broadening of the feature that indicates a combination of 
multiple peaks is not seen with the 13 C fraction. Thus 
the peak shift cannot be interpreted in terms of a combi- 
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nation of the components corresponding to 12 C and 13 C. 
The present results indicate that the vibration modes in 
the QCC are not very localized, but involve rather large 
molecular structures. The peak shift is related to the num- 
ber of carbon atoms involved in the vibration mode. The 
larger shift than the calculation suggests that more than 
two carbon atoms are associated with the 6.2 /im band of 
the QCC. The large shift in the 11.4/iin band also indi- 
cates that this feature originates from CH bonds associ- 
ated with large benzene structures if it is ascribed to a 
CH out-of-plane bending mode. 

The two kinds of the QCC have very different chem- 
ical compositions from each other. The spectrum of the 
dark-QCC shows a continuum component, which indi- 
cates the development of sp 2 carbon-carbon bonds. The 
broad bump around 8.2 /im and the strong continuum 
are the characteristics of carbon-rich particles. The filmy- 
and dark-QCC have absorption bands at similar wave- 
lengths, but the bands seen in the dark-QCC always peak 
at longer wavelengths than the filmy-QCC except for the 
3.3 fim band with large 13 C fractions (Fig.HaJ). The isotopic 
shift appears similar both in the filmy-QCC and dark- 
QCC except for the 3.3 /im band. The 3.3 /urn band in the 
dark-QCC is rather weak and thus further examination is 
needed to understand the different behavior of the 3.3 /im 
band between the filmy- and dark-QCC. 

4. Discussion 

The infrared spectra of the QCC show several similari- 
ties to the observed UIR bands, but they are not exactly 
the same. The peak wavelength of the 3.3 /im band of the 
f-QCC is in agreement with observations (Sakata et al. 
199dJ), while the exact peak wavelengths are not at the 
right positions for some other bands and there are dif- 
ferences in the band profiles. The 11.4/im band in the 
QCCs is obviously located at a longer wavelength than 
observed (11.23-11.25 /im). Some UIR bands, such as the 
6.2 /im and 11.2 /im bands, are not symmetric, but skewed 
to longer wavelengths, whereas most bands in the QCC 
seem to be more or less symmetric. Therefore a direct com- 
parison of the QCC spectra with observations is not very 
straightforward since the QCC is probably not the very 
material that exists in interstellar space. However we be- 
lieve that the QCC contains many of the structural units 
corresponding to the material in the interstellar medium 
because of the spectral similarities of the QCC and the 
UIR bands. In the following we examine whether or not 
observed variations of the UIR bands seen in post-AGB 
stars and PNe can be accounted for by the isotopic shifts 
observed in the QCC. We concentrate on the relative shifts 
in the band peaks of the QCC and the absolute peak wave- 
lengths and band profiles are not discussed. The results 
of the filmy-QCC are compared with the observations for 
clarity because the difference in the relative shifts between 
filmy- and dark-QCC is small except for the 3.3 /im band. 

Peeters et al. ( 20021) have investigated the 6-9 /im UIR 
bands in various objects and found that there are at least 



three classes, designated as A, B, and C, according to the 
peak wavelengths of the emission bands. Class A is a major 
class and includes a wide range of objects, such as H n re- 
gions, non-isolated young stellar objects (YSOs), reflection 
nebulae, and galaxies. They do not show appreciable vari- 
ations in the 6.2, 7.7, and 8.6 /im band features. Class B 
consists of isolated YSOs, post-AGB stars, and PNe. They 
exhibit quite a large variation in the peak wavelengths of 
the 6-9 /im features. The peaks are all shifted to longer 
wavelengths compared to class A objects. The observed 
ranges of the peak wavelengths are listed in Table ^ The 
apparent variation in the 7.7/im complex comes partly 
from the change in the relative strengths of the subcom- 
ponents. The peak shift in the 7.8 /im subcomponent is 
well observed, while that in the 7.6 /im subcomponent is 
less clear. Tabled indicates the observed range of the peak 
wavelength of the 7.8 /im subcomponent. Class C consists 
of only two objects (CRL2688 and IRAS 13146-6243). 
They show a very different spectrum from the other two 
classes. They have a band at 6.3 /im but do not show any 
features around 7.7 and 8.6 /im. Instead they show a broad 
band at 8.22 /im. Hence class C objects seem to have dif- 
ferent kinds of the band carriers from those in classes A 
and B. In the following discussion the peak wavelength of 
the 6.3 /im band in class C objects is not included in the 
observed variation since it will probably not be directly 
related to the isotopic effects. Although there are some 
internal variations present, the class in the 6-9 /im bands 
is well defined and almost all class B objects show the 
peaks shifted to longer wavelengths in all the 6.2, 7.7, and 
8.6 /im bands. 

Compared to the 6-9 /tm spectrum, the 3.3 and 
11.2 /tm UIR bands show small variations (Hony et al. 
2001), but the variations in the peak position do exist 
(e.g. Nagata et al. IT5%£1 Tokunaga et al. H^SSI Roche et 
al.EiSl). Tokunaga et al. Ijl991|) have suggested that there 
are two types of the 3.3 /im band from a small number of 
objects. Type 1 is a major class, in which the peak is lo- 
cated at 3.289 /im, while type 2 has the peak at 3.296 /im. 
They also suggest that type 2 objects have a narrower 
width than type 1. van Diedenhoven et al. I|2003J) have 
extended the investigation of the 3.3 and 11.2 /im UIR 
bands on the same sample as in Peeters et al. H2002fl . They 
found that there are also two classes in each band, desig- 
nated as A3.3, -B3.3, An. 2 , and -B11.2, respectively. Class 
A3. 3 corresponds to type 1 and -B3.3 to type 2 in Tokunaga 
et al. I)1991|l . A3.3 and An. 2 are major classes, and -B3.3 
and -B11.2 classes show the bands at 3.3/im and 11.2 /im 
shifted to longer wavelengths compared to A3. 3 and An.2, 
respectively. The observed ranges of the shifts in the 3.3 
and 11.2 /im bands are also shown in Tabled The shift 
in the 3.3 /tm band is quite small and that in the 11.2 /im 
band is modest. In general objects classified as class B in 
the 6-9 /im spectrum are also classified as B3.3 and -B11.2, 
but the correlation is not very good. Objects classified as 
class B in the 6-9 /im spectrum sometimes do not show de- 
tectable shifts in the 3.3 and 11.2/tm bands (e.g. BD +30 
3639) and vice versa (e.g. He2— 113). van Diedenhoven et 
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al. I|2003[) concluded that the appearance of the 3.3 and 
11.2 /im band variations is not tightly correlated with the 
UIR bands in 6-9 /im. 

The present results suggest that if 12 C/ 13 C < 10 the 
isotopic shifts are detectable particularly in the 6.2 and 
7.7 /im bands. A low 12 C/ 13 C ratio of this range is often 
suggested in post-AGB stars and PNe, and thus part of the 
peak shifts in these objects could be attributed to the iso- 
topic effects. The shift in the 11.2 /im of a medium degree 
can also be expected, while the shift in the 3.3 /im should 
be small. This shift pattern is qualitatively compatible 
with observations except for the 7.8 /im band. The 7.8 /im 
subcomponent shows a larger shift in its peak wavelength 
compared to the QCC sample. The 7.7 /im complex ex- 
hibit complicated profile variations and there may be other 
weak subcomponents present. The observed range of the 
3.3 /im, 6.2 /im, and 11.2 /im bands suggest that the ob- 
jects with the largest shifts should have 12 C/ 13 C ~ 2. 

We have searched in literature for the 12 C/ 13 C ratio 
of objects that had been observed by ISO/SWS. Tabled 
lists the objects and the SWS observing modes. SWS01 
was the full grating scan mode and the spectral resolution 
depends on the scan speed. SWS06 was the fixed-range 
grating scan mode and provided the highest spectral reso- 
lution in the grating mode of SWS (Leech et al. 120021). We 
used the highest spectral resolution data as much as possi- 
ble if there are several observations made for the same ob- 
ject. IRAS 23133+6050 is included in the list as a reference 
of the class A object. The Off-Line Processing data of ver- 
sion 10.1 (OLP 10.1) were obtained from the ISO Archival 
Data Center and reduced further by the Observers SWS 
Interactive Analysis Package (OSIA) version 3.0. 1 The 
peak wavelength for each band is estimated after the lo- 
cal continuum has been subtracted (see FigEJ). The peak 
position is slightly dependent on the assumed local con- 
tinuum. Strong Pf<5 (3.29699 /im) has been removed in the 
estimation of the peak wavelength of the 3.3 /im band for 
IRAS 23133+6050, NGC7027, and BD +30 3639. Tabled] 
lists the objects with the derived isotope ratio and their 
observed peak wavelengths. The peak wavelengths are in 
agreement with previous investigations within the uncer- 
tainty (Nagata et al. 113551 Tokunaga et al. ITMT1 Peeters 
et al. 12002(1 . The objects are listed in the order of the 
peak wavelength of the 6.2 /im band. Each band spectrum 
is normalized at the peak after the continuum has been 
subtracted and shown in Fig0 

Very few data of the carbon isotope abundance are 
available for post-AGB stars. Observations of 13 C isotopic 
lines of C2, CN, and CH + in several post-AGB stars re- 
sulted in negative detection of molecular features of 13 C 
except for one object (Bakker et al. IT39T1 Q]H§. They 
thus provided only lower limits for the 12 C/ 13 C ratio of 
> 10 — 20. The carbon isotope ratios in a few post-AGB 
stars and several PNe have been obtained mostly from ra- 

1 The OSIA is a joint development of the SWS consortium. 
Contributing institutes are SRON, MPE, KUL and the ESA 
Astrophysics Division. 



Table 2. Objects and the ISO/SWS observation mode 



Name 


Type 


TDT" 1 


Obs. mode 1 ' 


IRAS 23133+6050 c 


H n region 


56801906 


SWS01 (2) 


IRAS 21282+5050 


post-AGB 


15901777 


SWS01 (3) 


NGC 7027 


PN 


33800505 


SWS06 


BD +30 3639 


PN 


86500540 


SWS01 (3) 


IC5117 


PN 


36701824 


SWS01 (1) 


HR4049 


post-AGB 


17100101 


SWS01 (2) 


HD 44179 


post-AGB 


70201801 


SWS01 (4) 


CRL 2688 


post-AGB 


33800604 


SWS06 



a Target Dedicated Time (TDT) of the ISO observation. 

b Number in brackets indicates the scanning speed of the 

SWS01 mode. 

c Reference object of class A. 



dio observations of CO molecules. The intensity ratio of 
12 CO to 13 CO emission is thought to be a good measure 
for the isotope ratio if both lines are optically thin since 
the selective photodissociation of 13 CO should be coun- 
terbalanced by the charge exchange reaction (e.g. Mamon 
et al. I1988[) . However 12 CO lines often become optically 
thick, only lower limits of 12 C/ 13 C being estimated. It 
should also be noted that observations of CO probe the 
entire region of the circumstellar envelope. The post-AGB 
is a transient phase and a rapid change of the abundance 
is expected to occur in their envelopes. Thus the isotope 
ratio in the vicinity of the star might be different from 
that inferred from the CO intensity ratio. 

IRAS 21282+5050 is a post-AGB star and classified as 
class A in the 6-9 /im spectrum and ^3.3. The 3.3 /im band 
peaks around 3.288-3.290 /im (e.g. Nagata et al. [1988). 
Likkel et al. (1988) obtained that the intensity ratio of 
12 CO / 13 CO ( J=l-0) is about 100 and inferred 12 C/ 13 ~ 
200, taking account of the optical thickness effect of 12 CO. 
Based on the large gradient velocity model analysis, Balser 
et al. <|2002[1 obtained 12 C/ 13 C > 32 from observations 
of CO J=2-l and J=3-2 transitions. The observed peak 
wavelengths are compatible with the large 12 C/ 13 C ratio. 

NGC7027 is a bright young PN. It shows the class A 
6.2 /im band feature, while the 7.7 and 11.2 /im bands are 
slightly shifted to longer wavelengths relative to class A. 
Also the 7.8 /im component is stronger than the 7.6 /im 
component, indicating a signature of the class B object. 
The 3.3 /im band peak appears at the normal (class A) 
wavelength. This object has been observed by several 
groups in 13 CO molecular lines. Because of the optical 
depth effect of 12 CO the CO line intensity ratio should 
be taken as lower limits for the isotope abundance ratio 
of 12 C/ 13 C. They range from 65 (CO(J=1-0); Kahane et 
al. IT^2) to a recent value of 11 (CO(J=2-l); Josselin 
& Bachiller l2UfJ3"|) . suggesting that the 12 C/ 13 C ratio in 
the CO envelope of NGC7027 is probably not less than 
10. Except for the 7.7 /im band appearance the suggested 
range of 12 C/ 13 C is compatible with the class A classifi- 
cation. 
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Table 3. Peak wavelengths and carbon isotope ratios of the objects observed with ISO/SWS 



Object 


3.3 /xm 


Peak wavelength (/xm) 

6.2 (im 7.8 /xm 11.2 /an 


12 c/ 13 c 


Reference 6 


IRAS 23133+6050 


3.289 


6.219 


7.795 


11.218 






IRAS 21282+5050 


3.290 


6.219 


7.817 


11.245 


> 200, > 32 


1, 2 


NGC 7027 


3.292 


6.222 


7.830 


11.247 


> 31, > 65, > 25, > 11 


3, 4, 5, 6 


BD +30 3639 


3.292 


6.240 


7.850 


11.237 


> 4 


7 


IC5117 


3.295: 


6.25: 


7.82: 


11.28: 


23, 14 


6, 7 


HR4049 


3.296: c 


6.256 


7.880: 


11.285 


~ 1.5" 


8 


HD 44179 


3.296 


6.267 


7.861 


11.250 


> 22, > 12.5 


3, 9 


CRL 2688 


3.297 


6.290 






> 19, 32j^°, > 3, 20±*, 5 


3, 4, 5, 10, 11 



a Intensity ratio of CO2 bands (Cami & Yamamura 2001 ) 
b References 

1: Likkel, et al. fK%5l : 2. Balser et al. 1120021 : 3: Bakker et al. (fTW7t : 4: Kahane et al. (TIMSt : 5: Bachiller et al. flWl : 6: 
Josselin & Bachiller iTSUUSll ; 7: Palla et al. J2T5001 : 8: Cami & Yamamura 9: Greaves & Holland RWH : 10: Wannier & 

Sahai (tTWl : 11: Jaminet et al. RMH . 
c Doubly peaked. 




Wavelength (jim) Wavelength Qjun) 

Fig. 7. UIR band spectra of the objects with the measured 12 C/ 13 C ratio, a) the 3.3 /xm band, b) 6.2 /im band, c) 
7.8 //m band, and d) 11.2 /xm band. A: IRAS 23133+6050, B: IRAS 21282+5050, C: NGC 7027, D: BD +30 3639, E: 
IC5117, F: HR4049, G: HD 44179, and H: CRL 2688. The spectra are normalized after the local continua have been 
subtracted. Strong line emissions have been removed for clarity. 
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BD +30 3639 is a young PN and classified as class B 
in the 6-9 /im spectrum. Palla et al. f20U0| reported a 
tentative detection of 13 CO (J = 2 — I) and obtain a lower 
limit of 12 C/ 13 C as > 4. This object shows medium size 
band shifts, which can be accounted for by the isotopic 
shifts with 12 C/ 13 C ~ 5. The corresponding shift in the 
3.3 /im would be only ~ 0.03 /im, which is also compatible 
with the observed peak position. The 11.2 /im band shows 
a large tail in the longer wavelength side. 

Assuming that CO lines are optically thin, Palla et 
al. <|2DUD}| derived the 12 C/ 13 C of 14 for IC5117, while 
Josselin & Bachiller {(20031 indicated the ratio of 23 in this 
object. IC5117 was not included in the sample of Peeters 
et al. H2002fl . The SWS spectrum is noisy and the peak 
wavelengths cannot be determined accurately. However it 
clearly indicates the characteristics of the class B spec- 
trum. The 7.7 /im band complex is dominated by the 
7.8 /im component and the 3.3 and 6.2 /im bands peak at 
longer wavelengths than class A. The indicated 12 C/ 13 C 
ratio seems to be too small to account for the observed 
spectrum by the isotopic effects, although a qualitative 
comparison is difficult because of large uncertainties in 
the peak wavelengths. 

HR4049 is a very metal-poor post-AGB star and it is 
difficult to estimate its carbon isotope ratio because of 
very few metallic lines detected in its spectrum. Recently 
Cami & Yamamura (2001) have investigated circumstel- 
lar CO2 features in the infrared region and suggested that 
this star has a very peculiar oxygen isotope ratio, such 
that 17 and 18 are quite enhanced. They have detected 
13 C 16 2 , 16 13 C 18 0, and 16 13 C 17 C lines, and the in- 
tensity ratios of the 13 C to 12 C lines were all about 0.6- 
0.7. Although detailed analysis is required to derive a re- 
liable ratio, this implies that HR4049 is quite rich also in 
13 C. The positive correlation of 12 C/ 13 C with 16 0/ 17 
seen in several environments (Wannier & Sahai[I2H21 also 
suggests a low 12 C/ 13 C ratio. CO2 molecular emission 
originates from the region near the star and thus the de- 
rived ratio should correspond to the region in which the 
UIR bands are also emitted. The observed UIR bands in 
HR4049 are all shifted to longer wavelengths relative to 
normal class A positions, which can be accounted for by 
isotopic shifts with 12 C/ 13 C - 3. The small 12 C/ 13 C ra- 
tio suggested by the CO2 line emissions is compatible with 
the observed shifts. 

HD44179 is a well-studied post-AGB star and very 
bright in the infrared. Greaves & Holland ifTMTjl made ob- 
servations of CO (J=2— 1) and obtained the peak intensity 
ratio of 2.2 ± 0.2 for 12 CO/ 13 CO. The central part of the 
12 CO profile appears optically thick and from the wing of 
the profile they estimate 12 CO/ 13 CO> 12.5. The negative 
detection of 13 CH+ suggests 12 C/ 13 C > 22 (Bakker et al. 
I1997|l . The 6-9 /im spectrum of HD44179 belongs to class 
B and shows one of the largest shifts in the peak wave- 
lengths among the class. HD44179 is also the only one ob- 
ject that shows two subcomponents in the 6.2 /im clearly. 
The 3.3 /xm band of HD44179 is of type 2 and peaks at 
3.296 /im. The 11.2 /im band is also shifted to 11.25 /im. 



These shifts can be accounted for consistently by the iso- 
topic effects if 12 C/ 13 C is ~ 2. This required ratio is ob- 
viously smaller than the one indicated by CO and 13 CH + 
observations. A complicated torus-bipolar cone structure 
resulted from a transition from oxygen-rich to carbon-rich 
envelope has been suggested in HD44179 (Waters et al. 
IH?9"8l Men'shchikov et al. l2"0U2) . indicating that the car- 
bon isotope ratio may also spatially vary. In fact Kerr et al. 
( 199(2J) have indicated that the 3.3 /im band peak changes 
with positions within the nebula and that the longer wave- 
length peak appears near the star and the interface re- 
gion of the biconical structure. The 13 CO emission with 
a similar strength to 12 CO has also been detected for the 
4.6 /xm bands in the nebular region (Waters et al. Q998). 
Two components seen in the 6.2 /im band in HD44179 may 
then be ascribed to the band carriers with a small 13 C iso- 
tope abundance in the outer part and to those with a low 
12 C/ 13 C formed in the region close to the star. Unless 
there is an appreciable increase in the 13 C abundance in 
the vicinity of the star, however, the observed peak shifts 
in HD44179 cannot be accounted for by isotopic effects 
alone. High-spatial resolution observations in other UIR 
bands are interesting to investigate possible spatial varia- 
tions in their peaks (Song et al. 2003). 

CRL 2688 is a carbon-rich post-AGB star which shows 
a peculiar 6-9 /im spectrum (class C). The integrated 
intensity ratio of 12 CO to 13 CO (J=2-l) is about 3-4 
(Wannier & Sahai 113571 Bachiller et al. lTW7|l . The 12 CO 
line is optically thick and thus this gives a lower limit. 
Wannier & Sahai <|1987|1 derived the CO/ 13 CO ratio to 
be 20 from model analysis on CO (J=2-l) observations 
and Kahane et al. 1)1992(1 obtained the ratio of 32^° from 
observations of CS. Jaminet et al. (1993) have indicated, 
on the other hand, that the fast wind component has 
12 C/ 13 C of 5 from the wing profile of CO (J=3-2) emis- 
sions. Bakker et al. gave a lower limit of 19 based 
on the negative detection of circumstellar molecules with 
13 C. Goto et al. l|2HU2|l indicated that the 3.3 /im band pro- 
file does not show any spatial variations, suggesting that 
the extended component of the band emission is the scat- 
tered light of the emission from the central region. There 
is no evidence so far for the spatial variation in the band 
features within the nebula. Since this object may have a 
different type of the band carriers, the peak position of 
the 6.2 /im band (in fact it is located at 6.29 /im) should 
be discussed separately from other objects. The 8.2 /im 
band seen in the dark-QCC may correspond to the band 
at 8.22 /im detected in class C objects. It should be noted 
that the dark-QCC of the natural isotope abundance also 
has a band at 6.3 /im, shifted to a longer wavelength com- 
pared to the filmy-QCC. The 3.3 /im band is seen at the 
class B position and the 11.2 /im is very weak, if present. 
Differences in the type of the band carriers themselves 
may account for the observed shift in the 6.2 /im band as 
suggested by the difference between filmy- and dark-QCC 
spectra. 

As described above, carbon isotope ratios obtained 
from CO observations of post-AGB stars and PNe are 
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upper limits in most cases or show large scatter because 
of uncertainties in the observed intensities of faint 13 CO 
emissions (see Josselin & Bachiller 2003) and in the cor- 
rection for the optical depth effects of 12 CO lines. In ad- 
dition very few objects in the post-AGB phase with the 
UIR bands have measured 12 C/ 13 C data, including upper 
limits. Thus it is difficult at present to directly compare 
the present results with observations for individual ob- 
jects. The overall shift pattern of the UIR band peaks is 
qualitatively in agreement with the isotopic shift pattern 
except for the 7.7/im band complex. However, Tabled 
does not indicate a quantitatively clear correlation be- 
tween 12 C/ 13 C and the peak wavelengths and the sug- 
gested 13 C abundance in individual objects seems to be 
slightly too small to account for the observed shifts in 
terms of the isotopic effects. Observations also indicate 
that the correlation among the peak wavelengths of the 
6-9 /jm bands is good, but the peak shifts in the 6-9 /im 
bands are not tightly correlated with those in the 3.3 and 
11.2 /im bands. The UIR band emission is thought to come 
from the vicinity of the central star, while the CO emis- 
sion probes the entire region of the circumstellar envelope. 
Part of the lack of clear correlations and the poor quanti- 
tative agreement between the isotope ratio and the peak 
wavelengths may be attributed to possible spatial varia- 
tions in the isotope abundance within the object. If the 
3.3 and 11.2 /im bands come from regions different (either 
closer to or farther away from the source) from that of 
the 6-9 /urn bands, the poor correlation of the peak wave- 
lengths among the 3.3, 6-9, and 11.2/xm bands may also 
be attributed to the spatial variation in the isotope ratio 
with which each band carrier was formed in the envelope. 

Most PAHs have band peaks longer than 6.3 /zm for 
a C=C vibration contrary to the filmy-QCC and it has 
been suggested that the substitution of a carbon atom by 
a hetero-atom, such as nitrogen, oxygen, or silicon, will 
shift the band to shorter wavelengths as observed in class 
A objects (Peeters et al. 12002J). Then the observed vari- 
ations are interpreted in terms of the relative abundance 
of pure-carbon PAHs and substituted PAHs. The hetero- 
atom substitution, on the other hand, does not make sys- 
tematic shifts in other bands. The isotopic substitution 
of 13 C has similar effects on the peak wavelength for the 
6.2/im band but in the opposite direction. In addition it 
makes systematic shifts in the peak wavelengths of other 
bands and qualitatively accounts for the observed varia- 
tions correlated with the 6.2/im band. The peak wave- 
length of the C=C stretching vibration band decreases 
with the size of PAHs, while the increase of the molecular 
size does not make systematic shifts in the 7.7/im band 
(Peeters et al. I2()()2|l . Hetero-atom substitutions and size 
variations are likely to occur in interstellar medium and 
the observed variations could be a summation of various 
processes. The present results indicate that the isotopic ef- 
fects can contribute to the observed variations for objects 
with low 12 C/ 13 C ratios. It should also be emphasized 
that the peak shifts of the UIR bands are observed mostly 



in post-AGB stars and PNe, which are the objects where 
small 12 C/ 13 C ratios are expected. 

It has been suggested that deuterated PAHs (PADs) 
can be efficiently formed in dense molecular clouds be- 
cause the large difference in the zero-point energy between 
H and D containing species leads to preferential formation 
(fractionation) of deuterated species at low temperatures 
(e.g. Tielens HW7\ Sandford et al. l2lTOjl . For the carbon 
isotope case, the zero-point energy difference is not large 
and the temperature in the circumstellar envelope is much 
higher then the zero-point energy. Hence isotopic fraction- 
ation should not be significant. Instead a rather large 13 C 
abundance is expected to exist in certain environments 
and carbonaceous dust that is formed in such environ- 
ments should show isotopic shifts in its band features. 
Shifts in the band peaks of PADs from those of PAHs 
are quite large because of the large mass difference in H 
and D. Contrary to deuterated species, expected shifts for 
13 C carbonaceous material are small but should still be in 
the detectable range as shown by the present results. 

5. Summary 

We synthesized the QCC, a laboratory analogue of car- 
bonaceous dust, with various 13 C fractions and measured 
the isotopic shifts in the peak wavelengths of the infrared 
bands. They all shift to longer wavelengths approximately 
linearly with the 13 C fraction. The fact that no separate 
peaks originating from 13 C-QCC appear indicates that the 
vibration modes of infrared bands associated with carbon 
atoms in the QCC should not be very localized, but that 
they stem from rather large structures containing sev- 
eral carbon atoms. The shifts in the 6.2 and 7.8 /im are 
quite large (AA > 0.2 /im per 13 C fraction) and that in 
the 11.2 /im is modest (~ 0.16 /im), while the shift in the 
3.3 /im band is small (< 0.02/im). The small shift in the 
3.3 /im band is in agreement with a simple calculation of 
benzene molecules and those in the 6.2 and 11.2/xm seem 
to be larger, also suggesting that a larger number of carbon 
atoms associated with these bands than benzene. The iso- 
topic shifts obtained in the present experiment are qualita- 
tively in agreement with the shift pattern observed in the 
peak wavelengths of the UIR bands except for the 7.7 /im 
band complex. However the observed variations seem to 
be larger than those inferred from the 12 C/ 13 C ratio and 
the quantitative agreement between the carbon isotope 
ratio and the peak wavelengths is not good for individ- 
ual objects. This may be attributed partly to large uncer- 
tainties in the isotope ratios derived from observations, 
to possible spatial variations in the isotope abundance in 
the envelope, and to combinations with other effects, such 
as hetero-atom substitutions. The present results indicate 
that the isotopic shifts in the peak wavelengths of the UIR 
bands should be detectable in objects of low 12 C/ 13 C ra- 
tios and part of the observed variations in the band peaks 
can be attributed to the isotopic effects. 
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Appendix A: Calculation of isotopic peak shifts of 
benzene molecule 

In this appendix, we estimate isotopic peak shifts in the 
vibration modes of a benzene molecule to compare with 
the isotopic shifts in the corresponding vibration modes 
of the QCC. Benzene is the simplest PAH and the results 
should be used as a 0-th order reference in comparison 
with the isotopic shifts in the QCC. 

The GF matrix method is used to calculate isotopic 
shifts of a benzene molecule (Wilson et al. |T955 ) . The G 
matrix represents the molecular geometry and the F ma- 
trix consists of the force constants. The eigenvalue of the 
GF matrix is given by A = 2ir 2 v 2 , where v is the fre- 
quency of the vibration mode. In the calculation of the 
isotopic shift, only the values of the G matrix elements 
change due to the isotopic substitution and the F matrix 
elements remain unchanged if all carbon atoms are sub- 
stituted by 13 C. We use the values of the force constants 
of benzene given by Crawford & Miller ifTMfl lIMty . The 
vibrational modes of infrared active fundamentals in ben- 
zene consist of one A 2u (C-H out-of-plane bending) and 
three E\ n (three kinds of in-plane vibrations) modes. The 
single G matrix element of A 2n is written by 

G(A 2u ) = (m + Hc) <J 2 (A.l) 

and those of E\ n are given by 

G(E lu ) 

§Mc -^(2a + 3r) M c 

fT 2 MH + (<T 2 + \<>T + §T 2 V C J 

(A.2) 

where only the upper right elements are shown for the 
symmetric G matrix. The parameters a and r indicate 
the reciprocals of the C— H and C— C bond lengths, while 
/zh and /xc are the reciprocal masses of H and C, respec- 
tively. The bond angles are all set as 120°. From Ecis. ljA.lfl 
and (|A.2|I and the F matrix we calculate the eigenval- 
ues of the GF matrix and obtain the wavelength ratios 
of each vibrational mode between 12 C-benzene and 1 C- 
benzene. Then the wavelengths of 12 C-benzene are scaled 
to the measured values. The results are listed in TablelA~T1 
Because the wavelength of each mode is different from that 
in the QCC, we scale the wavelengths of 12 C-benzene to 
the bands measured in the 12 C-QCC and estimate the 
isotopic shifts (Table The A20 mode is scaled to the 
3.3 /im band and the A19 mode to the 6.2 /mi band. The 
band at 11.4 /xm in the QCC is assigned to a solo or duo 
C— H out-of-plane bending, which does not have a corre- 
sponding mode in benzene. We estimate the isotopic shift 



Table A.l. Wavelength of the vibration modes of 12 C- 
and 13 C-benzene 



mode 


i;l C-benzene 

(H 


ia C-benzene 
(/mi) 


A20 


3.245 


3.257 


A19 


6.812 


6.964 


Ais 


9.524 


9.709 


An 


14.90 


14.95 



by simply scaling the An mode to 11.4 /xm because it is 
also a C— H out-of-plane mode of benzene. 
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